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We discuss the algorithm of the cutting rules of calculating the imaginary part of physical ampli- 
tude and the optical theorem. We ameliorate the conventional cutting rules to make it suitable for 
actual calculation and give the right imaginary part of physical amplitude. The calculations of the 
imaginary parts of several Feynman diagrams show that the ameliorated cutting rules agrees with 
the conventional integral algorithm very well. The investigation helps us to find that the optical 
theorem has severe contradictions and problem, thus isn't right. Besides, the calculation of a phys- 
, ical amplitude show that the ameliorated cutting rules keeps the physical result gauge-parameter 

q . independent. 
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I. INTRODUCTION 



It's well known that the discontinuity of a Feynman amplitude comes from the singularity of the Feynman propagator 
I/"") ' 0,0. In an 1-Particle-Irreducible (1PI) Feynman diagram the discontinuity will contribute imaginary part to the 
\ Feynman amplitude. One of the way to calculate the discontinuity is the causal perturbative theory Q, another is 
called as cutting rules proposed by Cutkosky .lj. Here we mainly discuss the cutting rules. We note that although 
the cutting rules is well-defined, it hasn't given an integrative algorithm to practically calculate the discontinuity. So 
we want to investigate it further. 

On the other hand, the optical theorem has given a strong constraint on the imaginary part of physical amplitude. 
Whether they agree with each other hasn't been clearly investigated. So we will discuss their relationship and 
investigate the origin of the optical theorem. The arrangement of this paper is: firstly we discuss how to ameliorate 
the cutting rules and investigate its physical meaning; then we calculate the imaginary parts of some Feynman 
diagrams to see if the ameliorated cutting rules coincides with the conventional integral algorithm; in section IV we 
investigate the breaking of the optical theorem; in Sect.V we calculate a physical amplitude to see if the physical 
^H. result keeps gauge independent under the ameliorated cutting rules; lastly we give the conclusion. 
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II. AMELIORATE THE CONVENTIONAL CUTTING RULES 



;_i ' We usually encounter branch cut in Feynman amplitude calculation. Such branch cut will contribute imaginary 
9^ i part to Feynman amplitude in an 1PI Feynman diagram. Consider the Feynman propagator 

D f(* -V)= [ ( t\ 2 V ■ e- iP<X - y) > (!) 
J (27r) 4 p z — m z + is 

there exist two singularities po — ±((p 2 + m 2 ) 1 / 2 — ie) in the denominator. They separately represent the processes of 
the on-shell particle propagating from position y to position x and from position x to position y. The cutting rules is 
just the algorithm to calculate the contribution of the processes to Feynman amplitude. But it isn't suitable for actual 
calculation since it assumes that only one of the two singularities of Eq.(l) has contribution to Feynman amplitude, 
however, it doesn't tell which singularity is the case In the follows we investigate this problem from the imaginary 
part of the unstable particle's self energy. We firstly calculate the imaginary part of the unstable particle's self energy 
when keeping all of the contributions of the singularities of the Feynman propagators to Feynman amplitude. At 
one-loop level the result is 



Im(—i) 



d 4 k 1 1 (m 4 + mf + m| — 2m 2 m\ — 2m 2 m 2 — 2m 2 m 2 .) 1 / 2 



(27r) 4 k 2 — m\ + ie (k — p) 2 — m 2 + ie Wwrn 2 

x (0[mi — m — m2\ + 0[m — mi — m,2[ + B\m,2 — m — mi]) ,(2) 



where the external-line momentum p is on shell po = (p 2 + m 2 ) 1 / 2 and 9 is the Heaviside function. According to the 
Breit-Wigner formula the imaginary part of the unstable particle's one-loop self energy is proportional to its decay 
width: ImM{p — > p) = mT Q- So one can easily see that the result of Eq.(2) is wrong since only the second term of 
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the right-hand side of Eq.(2) satisfies the Breit-Wigner formula. We can draw a sketch to illustrate the origin of the 
three terms of the right-hand side of Eq.(2), 




m 2 ra 2 m 2 



where the arrow denotes the propagator with it is cut and the momentum q along it satisfies the positive on-shell 
condition qo — (q 2 + M 2 ) 1 / 2 (M is the mass of the propagator). We call such momentum propagating direction 
denoted by the arrow as positive- on- shell momentum propagating direction. The three cut self-energy diagrams 
separately in turn represent the origins of the three terms of the right-hand side of Eq.(2). Obviously only the second 
cut is acceptable by the Breit-Wigner formula, the others must be eliminated. Comparing the three cuts we find 
that the only difference between them is that the two positive- on- shell momentum propagating directions of the cut 
propagators of the second cut are inverse in the momentum loop, but that of the others are equi-directional in the 
momentum loop. We can use this point to constrain which cut is acceptable, i.e. the positive-on-shell momentum 
propagating directions of the cut propagators must be inverse in every cut momentum loop. In an 1PI Feynman 
diagram such constraint makes all of the positive-on-shell momentum propagating directions of the cut propagators 
are along a same direction, i.e. the momentum-energy propagating direction of the external-line particles (see Fig. 1-5 
for a rough knowledge). We note that such picture can be used to describe the on-shell virtual particles' propagating 
processes happening in the quantum field vacuum, if such processes exist. Obviously in this picture the on-shell virtual 
particles satisfy the energy conservation law. Contrarily, if the above constraint doesn't exist, these on-shell virtual 
particles will not keep the energy conserved because they will provide with energy reciprocally in the cut momentum 
loops thus the total energy of them can be greater than that of the incoming particles. 

Obviously there are at most two cuts in each momentum loop based on this constraint. On the other hand, 
according to the knowledge that cutting once in a momentum loop is equivalent to performing the conventional loop 
momentum integral, each cut momentum loop can only be cut twice. For the uncut propagator, since its singularity 
has no contribution to Feynman amplitude, it should be replaced by its Cauchy principle value in the loop momentum 
integral. Summing up all the discussions we obtain the following ameliorated cutting rules: 

1. Cut through the Feynman diagram in all possible ways such that the cut propagators can simultaneously be put 
on mass shell and the cut propagators' positive-on-shell momentum propagating directions are reverse in every 
cut momentum loop; keep only two cuts in each cut momentum loop. 

2. For each cut propagator with definite positive-on-shell momentum propagating direction, replace l/(p 2 — to 2 + 
ie) — > — 2ir i 0[po] 6(p 2 — to 2 ) (where the momentum p is along the positive-on-shell momentum propagating 
direction), for each uncut propagator, apply Cauchy principle value to it, then perform the loop integrals. 

3. Sum the contributions of all possible cuts. 

One can easily find that the above ameliorated cutting rules is different from the conventional one in some aspects. 
Contrarily to the conventional cutting rules, the ameliorated one constrains that only two cuts can exist in one 
momentum loop and implies that both of the two singularities of Eq.(l) can have contributions to Feynman amplitude. 
Besides, it definitely points out that the Cauchy principle value should be applied to the uncut propagator in the loop 
momentum integrals. We note that to the best of our knowledge there exist no similar discussions before. But then 
there exist some calculations of the imaginary parts of the Green functions by the causal perturbative theory which 
coincide with our results Q (see section III). 

The above ameliorated cutting rules is very abstract. So we give some illustrations in Fig. 1-5. The arrows in 
Fig. 1-5 represent the positive-on-shell momentum propagating directions of the cut propagators. In order to explain 
the cutting conditions we give two examples in Fig. 2. 




FIG. 1: Cuts of the one-loop three-point irreducible Feynman diagram. 
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FIG. 2: Cuts of the one-loop four-point irreducible Feynman diagram. 




FIG. 3: Cuts of the two-loop two-point irreducible Feynman diagram. 




FIG. 4: Cuts of the two-loop three-point irreducible Feynman diagram. 




FIG. 5: Cuts of the four-loop two-point irreducible Feynman diagram. 

Since each cut contributes an imaginary factor i to Feynman amplitude, all of the cut diagrams in Fig. 1-5 contribute 
imaginary parts to the Feynman amplitudes. This agrees with the conventional knowledge. On the other hand, one 
can easily find that many possible cuts in the conventional cutting rules have been eliminated by the ameliorated one. 
For example all of the propagators in Fig.l and 2 can be simultaneously cut in the conventional cutting rules Q, but 
this is forbidden in the ameliorated one. One can also see that both of the two singularities of many cut propagators 
have contributions to the Feynman amplitudes (e.g. the vertical propagator in Fig.l). Besides, some of the cuts in 
Fig. 1-5 have been used to calculate the imaginary parts of the Feynman amplitudes 0, Q . 

III. COMPARISON OF CUTTING RULES AND CONVENTIONAL INTEGRAL ALGORITHM 

It's well known that the imaginary part of Feynman amplitude can be directly calculated by the conventional 
integral algorithm, i.e. by the Feynman parametrization, wick rotation and dimensional regularization In this 
section we will give four examples to see whether the results obtained by the ameliorated cutting rules agree with the 
ones obtained by the conventional integral algorithm. 
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Firstly we calculate the imaginary part of the two-loop two-point irreducible Fcynman diagram 
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where p 2 — m 2 . Using the conventional integral algorithm we have 
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where e = 4 — D, 7 is the Eulcr constant, and 
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In Eq.(3) we have used the formula In a — ln|a| — inO[—a]. On the other hand, using the ameliorated cutting rules 
we also have 
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Through numerical calculations we find that Eq.(5) is equal to Eq.(3). The concrete results are shown in Fig. 6. 
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FIG. 6: Results of Eq.(3) and Eq.(5). 



Secondly we calculate the imaginary part of the one-loop three-point irreducible Feynman diagram 
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where p 2 = (pi + P2) 2 = u„ and p\ — p\ — m 2 . Using the conventional integral algorithm we have 
d 4 k 1 1 1 
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where the mathematical formula (P denotes the Cauchy principle value) 

— - — = P- =F i7r<5(a) 
a ± ie a 



Vmi € (0,^), 
Vmi6(^t), (6) 

Vmi e (Sjsoo), 
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has been used in the Feynman parameter integrals. On the other hand, using the ameliorated cutting rules we also 
have (see Fig.l) 
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Obviously Eq.(8) is equal to Eq.(6). Eq.(8) also coincides with Eq.(3.24) of Ref . Q 

Thirdly we calculate the imaginary part of the one-loop four-point irreducible Feynman diagram 
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where ro = qo = 2p. Using the conventional integral algorithm and Eq.(7) we have 

T . , r d 4 k 1 1 1 1 

Im{— i) 
1 



16tt 2 



im / d 



(2ir) 4 k 2 + ie (k — q) 2 — to 2 + ie k 2 — m 2 + ie (k + r) 2 — to 2 + ie 
1 ,i-x f i-x- y d , 
dy 



Jo 



(3p 2 x 2 + 3p 2 y 2 — 10p 2 x y — 3p 2 x — 3p 2 y — to 2 z + m 2 — ie) 2 



6 



1 



167T 2 TO 2 



■1 fl-X 

Im I dx dy 
o Jo 

1 



1 



3p 2 x 2 + 3p 2 y 2 — 10p 2 x y + (to 2 — 3p 2 )x + (m 2 — 3p 2 )y — ie 



3p 2 x 2 + 3p 2 y 2 — 10p 2 x y — 3p 2 x — 3p 2 y + to 2 — ie 

r-1 rl-X 

dx I dy [5{3p 2 x 2 + 3p 2 y 2 - 10p 2 xy + (to 2 - 3p 2 )x + (m 2 - 3p 2 )y) 



1 



167rm 2 



o 



o 



— 5(3p 2 x 2 + 3p 2 y 2 — 10p 2 x y — 3p 2 x — 3p 2 y + to 2 )] 

_ °\P - f M7I - P] ln 9 P 3 -to 2 p + 2to 2 v /4 P 2 -to 2 
9p 3 — to 2 p — 2m 2 -\/4p 2 — to 2 



+ 



128tttoV 

m ~\n\2m 



In 



9p 3 — m 2 p + 2m 2 yj Ap 2 — to 2 



128tttoV 9(9 p 3 _ m 2 p _ 2m VV - ™ 2 ) 

^ ~ ^ ln - m 2 p + 2m W - m 2 )(15p 2 + 4 Pv /9p 2 - 12to 2 - 4to 2 ) 



1287rmV 9(9 p 3 m 2 p - 2m 2 ^/Ap 2 - m 2 )(15p 2 - 4 Pv / V - 12m 2 - 4m 2 ) 
On the other hand, using the ameliorated cutting rules we also have (see Fig. 2) 
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Obviously Eq.(lO) is equal to Eq.(9). 

Lastly we calculate the imaginary part of a little more complex Feynman diagram 
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where p 2 — M 2 and the cut diagrams have been shown. Using the conventional integral algorithm and Eq.(7) we have 
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On the other hand, using the ameliorated cutting rules we also have 
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Through numerical calculations we find that Eq.(13) is equal to Eq.(ll), as shown in Fig. 7. We note that a similar 
discussion has been present in Ref.0]. 

All of the four examples show that the ameliorated cutting rules agrees with the conventional integral algorithm 
very well. There are also many similar examples which we don't list here. 
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FIG. 7: Results of Eq.(13) (plotted as dot) and Eq.(ll) (plotted as line). 



IV. BREAKING OF THE OPTICAL THEOREM 



It's well-known that the optical theorem has put a strong constraint on the imaginary part of physical amplitude. 
Docs it really right? In this section we will discuss this problem. 

We firstly simply review the origin of the optical theorem. It's well known that the S-matrix is unitary: S'S = 1. 
Inserting S = 1 + iT one has: 

-i(T -Tt) = T*T. (14) 

Take the matrix element of this equation between the two-particle states IP1P2 > and |kik 2 >. To evaluate the 
right-hand side, insert a complete set of intermediate states, 

<PlP2 <PiP2|T t |/></|7 1 |kik 2 > , (15) 

/ / 

where the sum runs over all possible sets / of intermediate states and dY[f represents the phase space integral of /. 
Now express the T-matrix elements as invariant matrix elements M. times 4-momentum-conserving delta functions. 
Eq.(14) then becomes 



M(hk 2 ^ Pl p 2 ) - M*{j>xp 2 ^k x k 2 ) =iJ2 [ dl[M*( PlP 2^ f)M(hk 2 ^ f). 

f J f 



(16) 



Setting ki = pi and applying the kinematic factors required by Eq.(16) to build a cross section, one obtains the 
standard form of the optical theorem, 

ImM(k 1 k 2 — > hk 2 ) = 2E cm p cm a to t(kik 2 -» anything) , (17) 

where E cm is the total center-of-mass energy and p cm is the momentum of either particle in the center-of-mass frame. 
Considering the 1 — > 1 process one obtains the well-known relationship, 

ImM(a^a) = m a T a , (18) 

where m a and T a is the mass and decay width of particle a. 

It seems that the optical theorem is a straightforward consequence of the unitarity of the S-matrix. But we find 
that there are some severe contradictions in the optical theorem. We firstly discuss a concrete physical amplitude. 
It is known that the imaginary part of physical amplitude comes from two cases: one is the coupling constant has 
imaginary part, the other is the Feynman amplitude has branch cut. At one-loop level a physical amplitude can be 
expressed as 

M(i^ f) = ^2g k (a k + ib k ), (19) 



where the sum runs over different interactions, a k is a real number, 6 fc is the imaginary part coming from the branch 
cut of Feynman amplitude, and all of the imaginary parts of the coupling constants are included in the coefficient 
gk- For convenience we define two concepts: the quasi-imaginary part and the quasi-real part of Feynman amplitude, 
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which separately represent the parts of Feynman amplitude coming from the branch cut and normal integral. For 
example in Eq.(19) the quasi-imaginary part of the physical amplitude is ^2 k gkbk and the quasi-real part of the 
physical amplitude is ~^2 k g k a k . Under the charge conjugation transformation only the imaginary part of the coupling 
constant changes into its contra-value, all of the other quantities in Eq.(19) keep unchanged. So we have from the 
CPT conservation law 



M*(/->i) = M*(i^f) 



k 



g k (a k -ib k ) 



From Eqs.(19) and (20) we have 

M(i^f)-M*(f^i) = 2z]T 5fe & fe = 2ifmM(i^f), 



(20) 



(21) 



where Im takes the quasi-imaginary part of the physical amplitude. Obviously this result clearly demonstrates the 
relationship between the optical theorem and the cutting rules. If the optical theorem is right, one will have for the 
physical amplitude of top quark decaying into charm quark and gauge boson Z 



ImM(t^cZ) = \Y,I dl[M*(cZ^ /). 



(22) 



On the other hand, one can calculate the quasi-imaginary part of A4(t — > cZ) from the branch cuts of the Feynman 
amplitude. The one-loop Feynman diagrams of t — > cZ have been shown in Fig. 8. We note that there is no need 
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FIG. 8: One-loop diagrams of the physical amplitude of t — > cZ. 



to introduce counterterm into the one-loop physical amplitude of t — > c Z after including the contribution of the 
amplitude of Fig. 8. From the branch cuts of Fig. 8 we find a contradiction in the optical theorem: Eq.(22) forbids 
the cuts representing the contribution of f j d\\ ^ Ai* {c — > f)A4(t — > Z /) in which the cuts of t — > c two-point 
diagrams in the seven and eight diagrams of Fig. 8 are included, but that cuts of t — > c two-point diagrams in the 
seven and eight diagrams of Fig. 8 are permitted by Eq.(18). Obviously the contribution of the cuts of t — > c two-point 
diagrams in the seven and eight diagrams of Fig. 8 to the physical amplitude isn't equal to zero. 

In fact there exists a more severe contradiction. From the unitarity of S-matrix we also have T^T = TT^ . Thus 
Eq.(15) can be changed into 



<p lP2 |T t T|k 1 k 2 >= J d I[ <PiP2|r|/X/|Tt|k 1 k 2 > 



(23) 
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Since each intermediate state / has its own unique parameters (e.g. mass), Eq.(23) being equal to Eq.(15) means all 
of the terms about / are separately equal, i.e. for arbitrary / 

< Pl p 2 |Tt|/> = <p lP2 |T|/>, </|T|k 1 k 2 > = </|Tt| klk2> . (24) 

Since / are arbitrary, we ultimately obtain 

M(hk 2 - PlP2 ) - M*( PlP2 -» kxk 2 ) = . (25) 

Obviously Eq.(25) contradicts Eq.(16). 

Another problem about the optical theorem is: if the charge conjugation is conserved, from the CPT conservation 
law the left-hand side of Eq.(16) is a pure imaginary, but the right-hand side of Eq.(16) may contain real part which 
comes from the branch cut of the physical amplitudes. Why are there such contradictions and problem? In fact the 
S-matrix is only a symbolistic operator, it doesn't have the algebraic algorithm that a true operator should have, 
e.g. the multiplication and addition. This is because the S-matrix represents the summation of the interactions of 
the physical world, so it is no possible to express it by a concrete, having certain algebraic algorithm's, operator. 
Therefore the derivation of Eqs. (14-16) isn't right. 



V. GAUGE INDEPENDENCE OF PHYSICAL RESULT UNDER AMELIORATED CUTTING RULES 

In order to investigate whether the ameliorated cutting rules is reasonable, we calculate the physical amplitude 
t — > cZ to see if the ameliorated cutting rules keeps the physical result gauge independent. 

In Fig. 8 we have shown the one-loop diagrams of the physical process t-> c2. Since there is no tree-level diagram, 
it is no need to introduce counterterm to one-loop level. We note that in the following calculations we have used the 
program packages FeynArts, FeynCalc, FormCalc and LoopTools pj. Our calculations have shown that the quasi- 
real part of A4(t — > cZ) is gauge-parameter independent. So we only need to consider the quasi-imaginary part of 
Ai(t — > cZ). Using the ameliorated cutting rules we obtain (see also Fig.l) 

I~mM{t -► cZ)\t = 0, (26) 

where subscript £ takes the gauge-dependent part of the quantity. From Eq.(26) one can easily see that all of the 
unphysical cuts (e.g. t — > G + di) in Fig. 8 have been cancelled out in the physical amplitude by the ameliorated cutting 
rules. This coincides with the conventional knowledge |8|]. Obviously such result makes the decay width of t — > cZ 
gauge independent to two-loop level. 

By the way, from the above discussions we know that the optical theorem gives a contradictory result about the 
quasi-imaginary part of M(t — > c Z). If keeping the cuts of t — > c two-point diagrams in the seven and eight diagrams 
of Fig. 8 and for the first six diagrams of Fig. 8 only keeping the first kind cut of Fig.l required by Eq.(22), we obtain 



t v \ , \- V 2l V 3 * e A (3 - As 2 w ) (x c 



W c 
2\V2 




x [x\ - 2{i w + x d .i)x c + (£ w - Xd.if) 0[m c - m dti - M w vCw] , (27) 



where is the left-handed helicity operator, V 2 i and Vs% are the CKM matrix elements [Of , e is electron charge, sw 
and cw is the sine and cosine of the weak mixing angle, Myy and £w is the mass and gauge parameter of gauge boson 
W, m c and m,d,i are the masses of charm quark and down-type i quark, x c = m^/M^ v and Xd.i — m di /M^ v . This 
result is gauge dependent and will lead to the decay width of t — > cZ gauge dependent. One can easily see that the 
gauge-dependent terms of Eq. (27) come from the unphysical cuts of t — > c two-point diagrams in the seven and eight 
diagrams of Fig. 8. 



VI. CONCLUSION 



In order to calculate the contribution of the singularity of Feynman propagator to Feynman amplitude we investigate 
the cutting rules and the optical theorem. Wc ameliorate the cutting rules in order to make it suitable for actual 
calculation and give the right result of the imaginary part of physical amplitude. The calculations of the imaginary 
parts of several Feynman diagrams show that the ameliorated cutting rules agrees with the conventional integral 
algorithm very well (see Fig. 6, 7 and Eq.(6,8,9,10)). On the other hand, through careful investigation we find that 
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the optical theorem has severe contradictions and problem thus isn't right. Besides, the calculation of the physical 
amplitude t — > cZ shows that the ameliorated cutting rules keeps the decay width of the physical process gauge 
independent to two-loop level (see Eq.(26)). 

In the viewpoint of the conventional cutting rules the branch cut of Feynman diagram contributes the imaginary 
part to physical amplitude. Our investigation further finds that although the Feynman diagram doesn't represent 
the real physical process, it does provide a physical-like picture to describe the on-shell virtual particles' propagating 
processes happening in the quantum field vacuum if such processes exist. The ameliorated cutting rules not only gives 
a feasible algorithm to calculate the contributions of the virtual processes to physical amplitude, but also can help us 
to discover the deep-seated physical meaning of the quantum field theory. 
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